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We experimentally determine the equation of state of a laser cooled gas. By employing the Lane-Emden 
formalism, widely used in astrophysics, we derive the equilibrium atomic profiles in large magneto optical 
traps where the thermodynamic effects are cast in a polytropic equation of state. The effects of multiple 
scattering of light are included, which results in a generalized Lane-Emden equation for the atomic profiles. 

A detailed experimental investigation reveals an excellent agreement with the model, with a two-fold signif¬ 
icance. In one hand, we can infer on the details of the equation of state of the system, from an ideal gas 
to a correlated phase due to an effective electrical charge for the atoms, which is accurately described by 
a microscopical description of the effective electrostatic interaction. On the other hand, we are able map 
the effects of multiple scattering onto directly controllable experimental variables, which paves the way to 
subsequent experimental investigations of this collective interaction. 
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The concept of equation of state, the relationship be¬ 
tween thermodynamic state variables such as pressure, 
P, temperature, T, and volume, V, has evolved beyond 
the original formulation of Clapeyron for an ideal gas, 
PV = NKbT. In particular, it has been argued that 
there is a universal form for the equation of state for 
solids [1]. Matter at nuclear density may also allow for 
a description in terms of an equation of state, with im¬ 
plications on astrophysical observations of neutron stars 
[2-4]. In cosmology, an equation of state is expressed in 
terms of the ratio of pressure P to energy density p [5], 
and can be measured from observations on supernova dis¬ 
tances [6] . Even Einstein’s equation of general relativity 
can be derived from the Gibb’s relation 5Q = TdS, with 
5Q being the energy flux and T the Unruh temperature, 
allowing for the interpretation of Einstein’s equation as 
an equation of state [7]. More recently, with the advent 
of laser cooling, attention has been given to the thermo¬ 
dynamics of degenerate Bose [8, 9] and Fermi gases [10]. 

In astrophysics, the equation-of-state plays a central 
role in the study of stellar structure, where hydrostatic 
equilibrium condition of a polytropic gas under the grav¬ 
itational field and thermodynamic pressure leads to the 
celebrated Lane-Emden equation [11, 12]. Based on this 
approach, Chandrasekhar derived the mass limit for a 
white dwarf [13]. Moreover, the Lane-Emden formalism 
has also been used to test alternative theories of gravity 
in stars beyond the standard models of stellar structure 
[14]. 

In this Letter, we extend the Lane-Emden formalism 
to experimentally determine the equation of state of a 
laser cooled gas, by directly measuring the atomic den¬ 
sity profiles of large magneto-optical traps. In our exper¬ 
iment, the hydrostatic equilibrium condition is provided 
(in first order) by the balance between the harmonic con¬ 
finement - the analog of the gravitational force - and the 
thermodynamic pressure, cast in the form of a polytropic 


equation of state. However, when the cooling lasers tuned 
close to the atomic resonance, multiple scattering of light 
occurs and an additional collective interaction appears, 
due to the exchange of scattered photons with nearby 
atoms [15, 16]. In this regime, the atoms experience a 
Coulomb-like long-range interaction [17], therefore allow¬ 
ing to regard the system as an effective one-component 
trapped plasma [18]. Thus, the condition of hydrostatic 
equilibrium results in a generalized Lane-Emden equa¬ 
tion, which encompasses the joint effects of harmonic con¬ 
finement, thermodynamics and radiation pressure due to 
multiple scattering of light. Here, we provide the experi¬ 
mental evidence of the multiple scattering of light in the 
equation of state of a laser cooled gas and extracting the 
polytropic exponent by fitting the density profiles to our 
theory. A microscopical description of the interaction in¬ 
duced by multiple scattering is introduced to explain the 
experimental deviations from the ideal gas. 

A fluid description of a laser-cooled gas confined in a 
magneto-optical trap (MOT) may be introduced with the 
usual continuity and Navier-Stokes equations [18-20] 

^ -k V • {nv) = 0, (1) 


dv 

at 


-I- V ■ 'Vv = — 


VP 

mn 


El 

m 


F 

m ’ 


( 2 ) 


where n and v represent the gas density and velocity 
field, respectively, and m is the atomic mass. The col¬ 
lective interaction due to multiple scattering of light is 
determined by a Poisson-like equation 


V • F, = Qn, (3) 

where Q = {aji — o'l)o'lIo/c represents the square of the 
effective charge of the atoms [17], lo is the total inten¬ 
sity of the beams and c is the speed of light. Here, an 
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and ul represent the emission and absorption cross sec¬ 
tions, respectively [21]. The term Ft encompasses the 
usual MOT cooling and restoring force as determined by 
Ft = —av — Kr, corresponding to a damped-harmonic 
oscillator. The trapping potential is then assumed to 
be harmonic, which is a reasonable assumption for the 
anti-Helmholtz coils configuration. The spring constant 
is approximately given by k = B/hk = k[5,Iq,Is) 

where represents the Bohr magneton, a = a((5, Iq, Is) 
the friction coefficient from Doppler cooling, 5 is the laser 
detuning and Is is the atomic saturation intensity. For 
magnetic fields gradients of the order ~ 10 G/cm, devia¬ 
tions from the harmonic case are obtained for very large 
traps, namely i? ^ 1 cm. In our experiments, we use 
R 2 mm, so we may neglect anharmonic effects in the 
confinement. The thermodynamic effects are cast in the 
form of a polytropic equation of state 


P{r) = C^n(r)\ (4) 


where 7 is the polytropic exponent, = P{0)/n(0)'^ is a 
constant determined by the thermodynamic properties of 
the system, and P{Q) and n(0) the pressure and density 
at the center of the cloud. Introducing the condition of 
hydrostatic equilibrium in Eqs. (1), (2) and (3) yields 
the generalized Lane-Emden equation [20] 


where n(r) = n{0)0{r) and r = with = 

a typical scale in the system and 

uiq = k/ui the trap frequency. The dimensionless param¬ 
eter n = Qn{0)/Smu;Q = Wp/SwQ, with LUp = y/Qn{0)/m 
the equivalent plasma frequency [19], is the of ratio mul¬ 
tiple scattering (plasma) to the trapping forces. The sta¬ 
bility of the solutions can therefore be directly related 
to D, with stable solutions existing for 0 < D < 1 , as 
confirmed both numerically as by linear stability analy¬ 
sis [20]. The two limits correspond to distinct physical 
relevant scenarios. In one hand, for smaller traps with 
10 ^ ^ 10 ® atoms and large laser detuning, the dynamics 
is determined by the thermal effects and multiple scat¬ 
tering is negligible, corresponding to the limit when the 
atoms have no effective charge, Q 0 , or equivalently, 
D —>■ 0. In this limit, the equilibrium density profiles, for 
a spherically symmetric cloud, are given by 
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The case 7=1 and Ci = kBT (isothermal gas) simply 
corresponds to the Maxwell-Boltzmann equilibrium 


n(r) = n(0)e“^(’')/'=^^ = n(0)e“’'"/'^", (7) 

with R = y^2kBT/mujQ the 1 /e radius of the cloud. Later 
we shall see that the isothermal case is the most relevant 


solution in this limit. On the other hand, for very large 
traps with N > 10® atoms and small detuning, jdj < T, 
with r being the linewidth of the transition, the process 
of multiple scattering dominates, D —^ 1 , and thermal ef¬ 
fects can be ignored. In this case, a qualitative analytical 
solution can be found setting 7 —>■ 0 , yielding a step-like 
profile 


n{r) = n{0)Q{r — R), (8) 

with n(0) = 3mujQ/Q, R = the radius of the 

cloud and 0(r — R) the Heaviside function. These two 
limiting cases - the temperature-limited and multiple¬ 
scattering regime, respectively - are well known and have 
been reported by the early experiments [15, 22], although 
no relation with an equation of state has been established 
so far. Our experiments provide a quantitative measure¬ 
ment of the intermediate regimes, both theoretically and 
experimentally, on the equation of state of the gas, and 
its dependence on the effective charge Q. We shall also 
mention that a third regime may be possible, namely 
the two-component regime [22]. In this case, there is a 
strong confinement near the center of the cloud, due to 
the influence of the magnetic field on the optical pump¬ 
ing between the Zeeman sublevels of the ground state, 
and a weaker confinement in the outer region due to the 
Zeeman shift of the various excited sublevels. However, 
in large traps, almost all the atoms occupy the weak con¬ 
finement volume, the trap dynamics is essentially the one 
presented here and the presence of polarization gradients 
in the laser fields no longer influences the behaviour the 
of system. The influence of this regime is thus safely ex¬ 
cluded both from our theoretical model and experimental 
analysis. 

We now perform a detailed investigation of the equi¬ 
librium atomic profiles in a cold trap. Our experimental 
apparatus consists of a MOT [23], where ®®Rb atoms 
are collected from a dilute vapour in a background pres¬ 
sure of ~ 10“® Torr. Six independent trapping (and 
cooling) laser beams cross the center of the trap with 
beam waist of w ^ 4 cm, power per beam P ~ 40 
mW and wavelength A ~ 780 nm. The beams are not 
retro-reflected, thus avoiding feedback instability mech¬ 
anisms [24]. The trapping laser operates on the D2 line 
of ®®Rb (F = 3 —>■ P' = 4), and is red-detuned by 6, 
which can be precisely controlled by a double passage 
through an acousto-optic modulator (AOM). The transi¬ 
tion linewidth is approximately V/2'k = !> MHz. A mag¬ 
netic field gradient (VP) created with a pair of water- 
cooled coils in an anti-Helmholtz configuration (zero field 
in the center of the trap) generates a spatially dependent 
Zeeman split of the energy levels, yielding the restoring 
force of the trap. An additional repump beam, operat¬ 
ing on the hyperfine levels F = 2 —>■ F' = 3 of the D2 
line repopulates the trapping transition. The repump de¬ 
tuning is set by searching for the maximum fluorescence 
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signal, corresponding the larger number of atoms in the 
trap. We thus obtain a cold cloud with T ^ 100 jiK 
and N ^ 10^ to IV ^ 10^° atoms, depending on the laser 
detuning 5. A CCD camera collects the fluorescence sig¬ 
nal, illuminating the cloud with far from resonance light 
((5 = —4r), to avoid multiple scattering during the imag¬ 
ing process. In this way we measure the atomic distri¬ 
bution of the trap, integrated along the line-of-sight of 
the camera. Two additional CCD cameras, positioned in 
orthogonal directions, allow us to monitor the shape of 
the cloud. By using half-wave plates we independently 
control the intensity of the six trapping beams to achieve 
a spherically symmetric atomic ensemble. 

For each experimental condition, determined by the 
laser detuning 5 and magnetic field gradient VB, we av¬ 
erage the experimental CCD profiles over 30 realizations. 
One-dimensional profiles are obtained by cutting through 
a direction crossing the center of the cloud, whose coor¬ 
dinates are defined as the “center-of-mass” of the two- 
dimensional image. Each experimental profile is fitted 
with the general solution of Eq. (5), numerically com¬ 
puted with a fourth order explicit Runge-Kutta method, 
and integrated along one arbitrary direction (our system 
is spherically symmetric). The agreement between the 
experimental data and the theoretical model is excellent, 
for the whole range of experimental parameters inves¬ 
tigated here - see Fig. (1). By deacreasing the laser 
detuning, i.e. approaching the resonance, we clearly ob¬ 
serve a transition from a Gaussian to a flattened den¬ 
sity profile, also known as water-bag profile (which cor¬ 
responds to a paraboloidal curve, when integrated along 
the line-of-sight). It corresponds to a crossover from the 
temperature-limited regime, with 7 ^ 1 and SI ^ 0 , to 
the multiple scattering regime, with SI —>■ I. 

In what follows, we determine how the physical quanti¬ 
ties of the model, the dimensionless plasma frequency SI 
and the polytropic exponent 7 , scale explicitly as a func¬ 
tion of the experimental parameters along the crossover. 
For that purpose, we take several averaged realizations 
of the 2D profile from which one-dimensional profiles are 
extracted. By numerically fitting the latter to the gen¬ 
eral solution of Eq. (5), we plot SI and 7 against S and 
Vi?. The error bars thus correspond to the statistical 
standard deviations of a large number of fitted profiles, 
as depicted in Fig. (2). Note that we are not attribut¬ 
ing any error to the detuning 6, which is stable in our 
set-up. On the other hand, the laser power is not con¬ 
stant, with a 5 ^ 10 % drift that we can not overcome. 
This may indeed induce some fluctuations in the fitted 
parameters, although its significance should be mitigated 
by using large statistics. The size of the statistical error 
bars do, in fact, reflect the experimental “jitter” associ¬ 
ated with not only the laser power fluctuations, but other 
parameters drifts. 

As expected, when working closer to resonance, the ef¬ 
fect of multiple scattering becomes more important, not 
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FIG. 1. (color online) Integrated atomic density profiles (in 
false color code), and a one-dimensional horizontal cut pass¬ 
ing through the center of the cloud (white points). The cor¬ 
responding fitting curves to numerically generated solutions 
of the Lane-Emden equation in (5) are also displayed (black 
dashed line). The magnetic field gradient varies from top to 
bottom as VB = 7.5, 10, and 12.5 G/cm. From left to right, 
we vary the detuning as <5 = —3.2, —2.8 and —2.4 F. The er¬ 
ror bars here correspond to standard deviations of the CGD 
pixel values, taken over 30 realizations. 


only because of the higher value of the cross sections an 
and ctl, but also because the number of atoms in the 
trap also grows, increasing the probability of a photon 
being reabsorbed before leaving the system. At the same 
time, larger deviations from the ideal gas (7 = 1 ) are also 
observed, as the gas starts to behave like a (weakly) cou¬ 
pled one-component plasma. Very close to resonances, 
|d| ~ 2F, we observe the onset of some mechanical in¬ 
stabilities in the trap, characterised by an oscillatory be¬ 
haviour of the fluorescence signal. This effect has been 
observed by other authors and reported in the literature 
as a self-sustained instability [25] , which are related with 
a competition between the confining force of the trap 
and the increasing repulsive interaction associated with 
multiple scattering. As our model relies on dynamically 
stable regimes, we excluded data taken for jdj < 2F. 

At this point, it is pertinent to query about the de¬ 
pendence of the equation of state, as determined by the 
the polytropic exponent 7 , on the effective charge of the 
atoms. To unravel this dependence, we develop a sim¬ 
ple microscopical theory for the interactions in the sys¬ 
tem. We start by explicitly computing the total energy, 
U = Ut + Uc, with Ut = ^NKbT the usual thermal 
energy and Uc the energy associated with the Coulomb 
interactions and responsible for the deviations from the 
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FIG. 2. (color online) In the left panel, we plot the normalized 
plasma frequency fl = i^p/Swo as a function of detuning S. 
The right panel depicts the polytropic exponent 7 , also as a 
function of detuning <5. The different curves correspond to 
VS = 7.5 G/cm (red line), VB = 10 G/cm (black line) and 
VB = 12.5 G/cm (blue line). 


ideal gas. To compute Uc we start with Poisson equa¬ 
tion, V^^(r) = —-^n(r), with q = ^/€oQ the effective 
charge of the atoms and n{r) = u(0) assumed to be con¬ 
stant throughout the cloud - remember the water-bag 
solution - which allows us to explicitly derive an ana¬ 
lytical solution. After determining the electrostatic en¬ 
ergy of the system, we can derive the corresponding pres¬ 
sure resulting from the exchange of scattered photons as 
Pq = and the total pressure of the gas 

reads P = Pq + with Pq = uKeT the ideal 

gas contribution. Correcting the pressure in the form 
of a polytropic equation of state as investigated above, 
P = C^rP, yields 


7=1-1- 


2/3e 

1+r 


with ^ 


1 / m 

15 \47rno / 


A detailed derivation of this result can be found in the 
Supplemental Material [26]. Note that, by assuming a 
constant density distribution, we are overestimating the 
correction of the polytropic exponent. For that reason, 
we make the substitution N —>■ Af®® = aN, with a ex¬ 
pected to be close to unit, a <1. In fact, allowing a to 
be a free fitting parameter yields a ~ 0.8 as expected. 
We finally obtain a universal form for the correction of 
the polytropic exponent, which agrees very well with our 
theory - see Fig. (3). 

In conclusion, we experimentally determined the equa¬ 
tion of state of a gas of cold atoms in large atomic traps 
by fitting the density profiles with the solution of a gen¬ 
eralized Lane-Emden equation describing the hydrostatic 
equilibrium of the gas. By explicitly evaluating the en¬ 
ergy associated with the effective electrostatic interac¬ 
tion, we were able to explain how the polytropic exponent 
depends of the mean-field potential of the atoms, which 
in its turn results from the exchange of scattered photons 
by the atoms. Our findings constitute a strong quantita¬ 
tive evidence of the fact that a laser-cooled gas can effec¬ 
tively simulate and behave like a weakly correlated one- 



FIG. 3. (color online) Polytropic exponent, 7 , as a function of 
the universal parameter for VB — 7.5 Glcm (red points), 
VB = 10 G/cm (black points) and VB = 12.5 G/cm (blue 
points). All points corresponding to different experimental 
conditions fall on the universal curve dehned by Eq. (9). The 
inset corresponds to the same experimental points but binned 
in the ^ parameter and equally spaced. 


component plasma. Recently, a Debye-Hiickel approach 
was introduced in the context of spin ice, where effective 
magnetic monopoles interact under a mutual Coulom- 
bic force [27], thus constituting another example where 
an effective plasmonic behaviour can be induced. The 
Lane-Emden formalism was, to the best of our knowl¬ 
edge, for the first time applied outside the context of 
astrophysics. The results obtained here pave the way to 
the subsequent experimental investigation of more exotic 
plasma-like processes in non degenerate cold gases. We 
have previously introduced the possibility of observing 
effects like phonon-lasing [28], classical rotons [29], plas- 
mon modes and Tonks-Dattner resonances, [18, 19], pho¬ 
ton bubbles [30], the dynamical Casiinir-effect [31] and 
twisted excitations carrying orbital angular momentum 
[32]. Multiple scattering should also play a role in the 
context of opto-mechanical instabilities in cold matter, 
recently proposed [33] and observed [34]. We conclude by 
referring to the close relation between the system investi¬ 
gated here and astrophysical processes involving trapped 
plasmas [35, 36]. Being able to achieve mimicking condi¬ 
tions in cold atoms laboratory experiments, with a great 
degree of control and tunability on the interactions, offers 
an ideal test bench to investigate astrophysical problems. 
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CPU/CQVibes, partially funded by EU FEDER, 
and from the EU FP7 projects LANDAUER (GA 
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Supplemental Material 

The purpose of this Supplemental Material is to guide 
the reader through the derivation of the correction of the 
polytropic exponent, 7. This parameter accounts for the 
deviation from an ideal gas (7 = 1) induced by multi¬ 
ple scattering of photons [15, 16]. The atoms behave as 
if they possess an effective electrical charge q = v^eoQ, 
with Q = {ur — cfl)ulIq/c [17], /q the total intensity 
of the beams and c is the speed of light. Here, an and 
gl represent the emission and absorption cross sections, 
respectively [21]. The induced collective interaction has 
been previously explored [18, 19]. 

Let us begin by determining the electrostatic poten¬ 
tial of the system, which is known to satisfy the Poisson 
equation 


V^(/)(r) =- qn{r). (10) 

eo 

We approximate the density distribution in the cloud by 
a water-bag profile - remember the multiple scattering 
regime discussed in the main text - corresponding to a 
constant density no spread over a radial extent of radius 
R, i.e. n (r) = n^Oir — R), with no = 'inuJ^jQ and 
/ \ 1/3 

R = i ) ■ This approximation allow us to keep the 

analysis tractable and derive an analytical correction for 
7. We shall then compute the solutions for the Poisson 
equation in two different regions. For outside the cloud, 
r > R, Eq. (10) reads, in spherical coordinates 

which admits solutions in the form ^ + 

B. Assuming that the potential vanishes at infinity, 
^(r—>-00) = 0 results in H = 0. Gauss theorem al¬ 
lows us to write A = with qx the total charge of 
the system, q^ = ^irquoR^. We then have, for r > R, 
^(r) = . Let now turn to the region inside the 

cloud, r < R, where the corresponding Poisson equation 
reads 





qnp 

eo 


( 12 ) 


In this case the solution are of the form (p (r) = 
Substituting in Eq. (12) results in A' = The 

integration constant B' is determined by the continuity of 
the potential (j){r) in the boundary of the two regions, i.e 
4'{R~) = Finally, we can write the electrostatic 

potential inside the cloud as 




6eo ^ 3eo 


(13) 


The next step is the evaluation of the effective electro¬ 
static energy, determined by 

t^C = ^ / (j){r)qn{r)dV (14) 

Introducing again the water-bag density profile yields 
the results Uc = ^ttQtiqR^ or, equivalently, Uc = 

i ypa terms of the volume and the number 

of particles in the system, which will be useful in the 
next steps. We now wish to evaluate the pressure in 
the cloud, which encompasses the contributions from the 
ideal gas part and the effective electrostatic interaction, 
P = Pp -\- Pc, with Pp = kBTn = ksTN/V and Pc 
determined by 


Pc = - 



(15) 


since the electrostatic energy doesn’t dependent on the 
temperature. We then have Pc = -^ 
or, equivalently, Pc = -^QripR'^. The total pressure in 
the system is given by 


ksTN 1 / 3 QN^ 

V F4/3 


(16) 


We now wish to establish an equivalence between the 
former equation of state and a polytropic-like one, in the 
form P = CjtP, as in the Lane-Emden derivation. With 
that in mind, we can write 


CjUp = ksTnp -\- 


Qn^R^ 


15 


or, equivalently, dividing by ksTrip 


C-i . QnpPi^ 


\hkBT 


(17) 


(18) 


where we defined e = 7 — 1. Note that we can rewrite this 
last expression in terms of the parameter of the model 
introduced earlier, namely the effective plasma frequency 
H scaling factor = p^^np. Simple 

mathematical manipulation of Eq.(18) finally yields 


7=1 + 


2/3^ 
C + 1 


(19) 


with ^ an adimensional universal parameter defined as 
^ where we use the total number of 

atoms as N = ^nupR^. 
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